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Abstract 

We carry out an exploratory study of the isospin one ao(980) and the isospin one-half k scalar 
mesons using Nf = 2 + 1 + 1 Wilson twisted mass fermions at one lattice spacing. The valence 
strange quark is included as an Osterwalder-Seiler fermion with mass tuned so that the kaon mass 
matches the corresponding mass in the unitary Nf = 2+1 + 1 theory. We investigate the internal 
structure of these mesons by using a basis of four-quark interpolating fields. We construct 
diquark-diquark and molecular- type interpolating fields and analyse the resulting correlation 
matrices keeping only connected contributions. For both channels, the low-lying spectrum is 
found to be consistent with two-particle scattering states. Therefore, our analysis shows no 
evidence for an additional state that can be interpreted as either a tetraquark or a tightly- 
bound molecular state. 



1 Introduction 



The Naive Quark Model (NQM) is - despite its simplicity - surprisingly successful in qualita- 
tively describing the experimentally observed meson and baryon spectrum. This success has led 
us to think of mesons and baryons as qq and qqq bound states, respectively. In particular, no 
mesonic state incompatible with the quantum numbers of a qq system has been confirmed, yet. 
However, there are a few exceptions [TJ [2] which cannot be described in the NQM. One promi- 
nent example is the Roper resonance, another not less prominent one is the presence of too many 
scalar states (i.e with quantum numbers J PC = ++ ) with mass below 2 GeV as compared to 
the expectation from the NQM. These scalars are the /o(600) or a, /q(980), /o(1370), /o(1500) 
and / (1710) with isospin 0, the i^(800) or re and i^(1430) with isospin 1/2, and the a (980) 
and ao(1450) with isospin 1. 

This excess of states compared to the NQM expectation suggests that the picture of mesons as 
qq bound states is too simplistic and it has to be complemented by other quarkonic and gluonic 
structures. Consequently, it has been speculated that some of these particles are tetraquark 
states, i.e. bound states of two quarks and two antiquarks, or predominantly gluonic in na- 
ture. For example, according to one favoured interpretation [2], the states /o (1370), /o (1710), 
Kq (1430) and ao(1450) might indeed have dominant qq components, as expected in the NQM but 
the state /o(1500) might be, predominantly, the lightest (0 ++ ) glueball [3], and the lightest of the 
scalar states might constitute a nonet with a dominant tetraquark contribution [U [H [8] . 
While such an interpretation is adopted by other authors, e.g. [9], there are also different sce- 
narios discussed in the literature, as for example in [10] 111]. 

Experimentally, many of the aforementioned scalar resonances are difficult to resolve as they 
have large decay widths and several decay channels that sometimes open up only within a short 
energy interval. The question whether a physical state is dominated by a qq, a tetraquark, 
a glueball or other hybrid wavefunctions is then typically investigated through the analysis of 
its production and decay modes. These are directly accessible in experiments and can be often 
measured rather accurately. This justifies the high experimental activity |11|I12] in investigating 
the composition of these states. It is thus crucial to develop a deeper theoretical understanding 
for the internal structure of these states. 

A theoretical understanding from first principles requires a non-perturbative method. Since 
Quantum Chromodynamics (QCD) is the theory of strongly coupled quarks and gluons, such 
a non-perturbative method is provided by lattice QCD. But investigating the states in lattice 
QCD is also a challenging endeavour: the distinction between scattering states, resonances and 
bound states is subtle on a Euclidean lattice with finite spacetime volume. In fact, there is no 
continuum spectrum in a finite spatial volume and the Hamiltonian has only discrete eigenvalues. 
In order to disentangle these different physical phenomena it is necessary to study the volume 
dependence of the discrete eigenvalues of the Hamiltonian [131 EH E3 HE]. In particular, the 
coefficients of the large volume expansion of the discrete eigenvalues are related to the phase 
shifts of the scattering process. Moreover, as the volume increases, the eventual resonances 
produce " avoided level crossings" of eigenvalues 0. 

This method requires the extraction of more than the ground state in a channel with given 
quantum numbers. These excited states are increasingly difficult to extract with sufficient 

Note, however, that also the threshold may display the same phenomenon |17| . 
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precision, even though the field has recently seen tremendous progress in the methodology. 
An additional complication is the appearance of fermionic disconnected contributions, which 
are notoriously noisy. Therefore, the available lattice results on scalar mesons and possibly 
existing tetraquark states are still limited (cf. e.g. [IH1 H3 [201 [2D E21 [23l EH [25]). Certainly 
more and in particular independent investigations are needed to gain a better understanding of 
these scalar states. 

In this paper we perform an exploratory study of the ao(980) and the k using Wilson twisted 
mass fermions. It is the first study of this kind with Nf = 2 + 1 + 1 dynamical quark flavours, 
using gauge configurations provided by the European Twisted Mass (ETM) collaboration [26 1 127 1 
[28| I29|. Note that in particular the dynamical strange quark might be important for studying 
scalar resonances. In this exploratory investigation we address the question, whether or not 
these states could be consistent with a tetraquark or a molecule interpretation. We focus on 
the precise computation of correlation functions of operators with quantum numbers of the 
ao (980) and k mesons using four quark interpolating fields ignoring fermionic disconnected 
contributions. The latter implies that there is no mixing among four quark, two-quark and 
gluonic states. Apart from obvious technical advantages and having a testbed of our method, 
there is another important reason for working in this approximation: in [20] bound states close 
to threshold in the 1 = and the 1 = 1/2 channels have been found in the same approximation. 
These bound states, found in addition to the expected scattering states, were interpreted as a 
possible indication for a tetraquark nature of the corresponding states. In our study, performed 
with a similar operator basis, but a different lattice discretisation, we do not observe such a 
bound state in the 1 = 1/2 channel. Moreover, we also do not observe it in the 1 = 1 channel, 
which was not considered in [20] . Note that parts of this work have recently been presented in 
a conference proceeding [30] . 

The reason for focussing on the ao and the k are the following: the ao has isospin 1 = 1, i.e. 
when choosing I z = ±1 only a single disconnected contribution is ignored. The k meson, on the 
other hand, mixes only with the K + tt channel. 

The paper is organised as follows: in section [2] we introduce the lattice formulation followed by 
a discussion of the operator basis in section O the results of our study are discussed in section [4] 
and we conclude in the last section. 

2 Lattice setup 
2.1 Lattice actions 

This work is based on gauge link configurations generated by the ETM collaboration |26|, [2"T], 
[28| [29] with the Iwasaki gauge action [31] and Nf = 2 + 1 + 1 flavours of twisted mass quarks. 

The light degenerate (u, d) quark doublet is described by the standard Wilson twisted mass 
action 1321 . 



while for the heavy (c, s) sea quark doublet the twisted mass formulation for non-degenerate 






X 
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quarks of [33] has been used, 

S h c, V ylx ih) ,X ih \U] = a 4 Y / X (h) (x)(D w (7n ) + i^ 5 T 1 + T 3f i 5 )x {h \x). 



In both cases Z?w denotes the standard Wilson Dirac operator, 

Av(mo) = ^(7m(v m + V*) -aV^V^) +m , (3) 

while x = (x iX ) an d X = (x^ c \x^) are the quark fields in the so-called twisted basis. 
For reasons explained in |34] the same value of the standard quark mass parameter tuq has been 
used in both sectors. 

When tuning the theory to maximal twist, automatic 0(a) improvement for physical quantities 
applies |33l 135] . This tuning has been done by adjusting itiq such that the PCAC quark mass 
in the light quark sector vanishes (cf. [28] for details). 

At maximal twist in a massless quark renormalisation scheme the renormalised quark masses 
are related to the bare parameters \i a and by 

mf = Zp 1 ^ - ^V-SJ > m c = Zp 1 (a <F + ^Vsj (4) 

[33j . where Zp and Zg are the renormalisation constants of the non-singlet pseudoscalar and 
scalar densities. In our simulations the values of \i a and fj-s have been adjusted by requiring that 
the resulting lattice kaon and D meson masses approximately assume their physical values \28\ 

m\ EH- 

For the computation of observables we use a twisted mass discretisation for valence s quarks, 
which is different from the sea s quarks ([2]). It is given by ([I]) with x^ ~ ► X^ = (x^ s \x^ s ^) 
and hi — > jjL s . We do this, to avoid the problem of mixing between s and c quarks, which 
is discussed in detail in [361 EZ]- Note that there are two possibilities to realize a valence s 
quark, x^ s an d X^ s \ which differ in the sign of the twisted mass term, ±i^t s 75. Strategies 
and consequences of choosing s + or s~ are discussed in detail in sections 13.1.21 and 13.2.21 The 
bare strange quark mass fis has been chosen such that kaon masses computed within this mixed 
action setup with flavour structure s + d and s~u (which are degenerate and known to have less 
discretisation errors than their s + u and s~d counterparts |40t |4"H I42j ) agree with kaon masses 
computed in the unitary setup [36] [37] , i.e. using (J2]) also for valence s quarks. 

In this work we consider six gauge link ensembles with simulation parameters given in Table [TJ 
They differ in the space-time volume (L/a) 3 x T/a and in the light u/d quark mass [ii- The 
lattice spacing a ~ 0.086 fm is the same for all ensembles. More details regarding these ensembles 
can be found in |28|. 



The discussion of meson and four-quark creation operators (cf. section [3|) and their quantum 
numbers is more convenient with quark fields in the "physical basis", (u,d) and (s + ,s~). This 
physical basis is related to the "twisted basis" 

(x {u) ,X {d) ) and (x {s+) ,X {s ] ) introduced in (pQ) 

and ([2]) according to 

( U \ - Jw**>/*( X H ^ ( S+ \ - jrwwpf X (S+) \ 

\d) - e U (d) ) ' \s- ) - e U ( * _) r ( ) 
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Ensemble /3 


(L/af x T/a 


/'■/ 


/At 


/'a 


/'■« 


OL 

(fm) 


(MeV) 


# of 
configs 


A30.32 1.90 


32 3 x 64 


0.0030 


0.150 


0.190 


0.02280 


0.086 


284 


672 


A40.32 


32 3 x 64 


0.0040 






0.02322 




324 


200 


A40.24 


24 3 x 48 


0.0040 






0.02300 




332 


1259 


A40.20 


20 3 x 48 


0.0040 






0.02308 




341 


500 


A50.32 


32 3 x 64 


0.0050 






0.02336 




362 


431 


A80.24 


24 3 x 48 


0.0080 






0.02328 




455 


1225 



Table 1: Gauge link ensembles considered in this paper. The notation follows [28l 



where oj is the twist angle, which we have tuned to maximal twist, i.e. oj = ir/2. 

When computing temporal correlation functions (Ot(t 2 )Ojfc(ti)), where Oj and O}. are e.g. meson 
or four-quark creation operators, we only consider quark propagators connecting time t\ and 
t2, but ignore propagation of quarks within the same timeslice, e.g. from t\ to t\. For mesons 
this amounts to neglecting so-called disconnected diagrams. For four-quark operators e.g. of 
tetraquark or two-meson type both singly and doubly disconnected contributions (cf. (b) and (c) 
of Figure [TJ are omitted. Consequences of not considering disconnected diagrams are discussed 
in the following sections. 




(a) Connected contribution. (b) Singly disconnected contribution. 




(c) Doubly disconnected contribution. 



Figure 1: Relevant contributions to a four-quark operator two point function. 

Finally it should be mentioned that at finite lattice spacing isospin and parity are not good quan- 
tum numbers in twisted mass lattice QCD. These symmetries are broken by terms proportional 
to 0(a). Consequently, when doing spectroscopy, one has to take mixing with states of opposite 
parity and different isospin into account. Of course, in the continuum limit these symmetries 
are restored and QCD is recovered. Mixing in the context of ao(980) and k is discussed in detail 
in sections 13.1.11 and 13.2.11 
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2.2 The pseudoscalar meson spectrum 

Since a pair of pseudoscalar mesons is rather light and can have the same quantum numbers as 
the scalar mesons ao(980) and k, these are the most relevant scattering states to consider. Their 
masses are approximately equal to the sum of the two masses of the corresponding individual 
mesons. Therefore, a precise and comprehensive knowledge of the meson spectrum is important 
for our analysis. 

As we ignore disconnected contributions, the r\ = uu + dd meson and the neutral tt = uu — dd 
become degenerate and there is an r//ry'-like meson with valence quark structure ss, but no light 
uu or dd valence quarks, which we denote by r] s . 

Another particularity stems from the valence action used for the strange quarks discussed above. 
The kaon and the rj s can be constructed using s + and/or s~ strange quarks, resulting in different 
values for the meson masses at finite value of the lattice spacing. Similarly, in Wilson twisted 
mass lattice QCD the neutral (connected-only) and charged pion mass values differ. 

All meson masses relevant for our investigation are collected in Table [2J 



ensemble 


m Tr{ud,du) 


m ir(uu-dd) 


m K(s+d,s-u) 


m K(s+u,s-d) 


m Vs (s+s-) 


m ri a (s+ s+ ,s~s~) 


A30.32 


284(1) 


494(6) 


576(7) 


704(2) 




876(1) 


A40.20 


341(2) 




599(3) 




774(2) 




A40.24 


332(1) 


530(7) 


593(1) 


723(2) 




882(1) 


A40.32 


324(7) 




588(5) 




779(2) 




A50.32 


362(7) 




601(9) 




783(2) 




A80.24 


455(1) 


625(3) 


635(1) 


753(1) 




885(1) 



Table 2: The pseudoscalar meson spectrum in MeV with disconnected diagrams neglected. 
Omitted mass values are not needed in the context of the tetraquark study presented in this 
paper. 



3 Creation operators and analysis details 

In this study we exclusively consider creation operators with four quarks (two quarks and two 
antiquarks). The structure of our four-quark operators is oriented at phenomenological expec- 
tations and ranges from four-quark bound states (molecules formed by two mesons and bound 
diquark-antidiquark pairs) to two essentially non- interacting mesons (two-particle operators). 

Of course, standard quark-antiquark operators, e.g. du for oo(980) and su for k, would also be 
of interest. However, since we neglect disconnected diagrams (cf. section [2~Tj) . such two-quark 
operators do not generate overlap to trial states created by four-quark operators. Consequently, 
in our setup four-quark operators and quark-antiquark operators probe different sectors, which 
is, why we do not consider the latter in the following. In a subsequent improved study we 
plan to include disconnected diagrams and to combine two- and four-quark operators in a single 
correlation matrix. 

Information regarding the used four-quark operators is summarised in Table [3j The operators 
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will be discussed in more detail below. 



ensemble 


L|U.ajl xv 


CP 1 1 ffp 


type 


Dirac structure 


smearing 


smearing 


a (980) 




A30.32 






KK molecule 


75, 7m > 7m 75 




A40.24 


Gaussian 


APE 


i] s ir molecule 
Kir molecule 


75 


— 

75, 7m, 7m 75 


A 80 94 






Hirni a ylr 

vJ.HJU.CH iv 




^75,^ 








molecule 


75 






no 


no 


diquark 
i-T + K two-particle 


C75 

75 


- 


A40.20 






rj s + 7r two-particle 


75 








molecule 


75,7/. 






Gaussian 


APE 


?/ s 7r molecule 
diquark 


75 










molecule 


75,7/. 




A40.32 


Gaussian 


APE 


?? s 7r molecule 
diquark 


75 




A50.32 


Gaussian 


APE 


KK molecule 
diquark 


75,7m 

C 7 5,C 





Table 3: Four-quark creation operators. 



3.1 Creation operators, a (980) sector (quantum numbers I(J P ) = 1(0 + )) 
The expected low-lying spectrum in the ao(980) sector (« 1000 MeV) is the following: 

• A two-particle 77 + ir and a two-particle 7/ + ir state. 

— In nature: 

* Mass 771(77 + vr) » 548 MeV + 140 MeV = 688 MeV [12]. 

* Mass 771(77' + tt) » 958 MeV + 140 MeV = 1098 MeV p2]. 

— In our lattice setup: 

* Due to neglect of disconnected diagrams 77 has flavour structure uu + dd and is 
degenerate with the neutral pion (cf. section [2T2]) : the 77 + tt state is orthogonal 
to any trial state obtained by using an operator containing s quarks, i.e. can be 
ignored in the following. 

* Due to neglect of disconnected diagrams 77' becomes 77^ (cf. section [272]) : masses 
772(773 + 7r ) ~ m ( 7 ls) + Tn(jr) depend on the gauge link ensemble and can be read 
off from Tabled 

• A two-particle K + K state. 

— In nature: mass m(K + K) 2 x 496 MeV = 992 MeV [12]. 
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— In our lattice setup: masses m(K + K) ~ 2m(K) depend on the gauge link ensemble 
and can be read off from Table [2j 



• Possibly a bound ao(980) state (might be of quark-antiquark, molecule or diquark-antidi- 
quark type), mass m(a (980)) = 980 ± 20MeV [E]. 



To be able to resolve these low-lying states, we consider the following operators: 



• Molecule type operators: 

^980)" = E ( S "( x )75^(x)) (d(xhs*(x)) (6) 

X 

= E( s "( x )^^(x))(^(x)75n(x)); (7) 

x 

since pseudoscalar mesons (mesons with spin structure 75) are the lightest mesons in a 
given flavour sector, one expects possible molecular bound states to be of pseudoscalar- 
pseudoscalar type. We also consider molecule type operators with 75 replaced by 7^ and by 
7j75. These operators enlarge our correlation matrices and allow us to study also excited 
states, in particular two-particle states with relative momentum (cf. section H|). 

• Diquark type operator: 



<Om = E( £a6C ^( x )^ 5 ci c ' T (x))( e ^n d ^(x)C775 S e (x)); (8) 

X 

diquarks with spin structure 75 are known to be the lightest [TJ [38], [39], which is, why we 
use 75 in this operator. We also consider diquark type operators with 75 replaced by 1. 
As before, the main reason is to enlarge our correlation matrices allowing us to study also 
excited states. 

Two-particle type operators: 

<( 9 for particle = (E^)75-(x))(E rf »75«(y)) 0) 

\ x y 

CoT^ 6 = (E s "( x )^ s ( x ))(E d »^-(y)); (io) 



these operators resemble states with two non-interacting mesons and, therefore, should be 
particularly suited to resolve two-particle K + K and i] s + it states. Note that terms with 
x = y in ([9]) and (|10p also appear in the molecule operators © and (J7J), which is, why 
two-particle K + K and rj s + tt states can also be resolved, even though only molecule 
and diquark operators are used. However, the generated overlap to two-particle states 
is significantly larger, when two-particle operators are applied, which in turn results in a 
signal of better statistical quality (cf. the numerical results in section B~T]) . 
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3.1.1 Mixing due to twisted mass symmetry breaking 

In twisted mass lattice QCD both parity P and isospin I are broken by 0(a). Consequently, 
one has to take mixing with states of opposite parity and different isospin into account. When 
there is mixing with rather light states (lighter than those one is interested in), problems arise: 
correlators are slightly contaminated by weakly decaying exponentials, which become dominant 
at large temporal separations, at which masses are usually determined. 

Before we discuss mixing due to twisted mass symmetry breaking, it is important, to understand 
the effects arising by neglecting disconnected diagrams in more detail. In such a setup the 
valence quark flavour structure is conserved, i.e. quark-antiquark pairs can neither be created 
nor annihilated. For ao(980) this implies that any state that mixes with the trial states created 
by our operators (cf. ([6]) to (jlOp ) must have valence quark flavour structure udss. 

This observation is particularly important, when discussing parity mixing, because at first glance 
there seem to be states of opposite (negative) parity, which are light, namely pions I{J P ) = 
1(0~). However, since pions only have a u and a d valence quark, but no ss pair, they are 
orthogonal to any state we probe with our four-quark operators. On the other hand udss four- 
quark states with negative parity are expected to be rather heavy, e.g. could be a pseudoscalar 
meson and a scalar meson like K + k. 

Since the z-component of isospin I z is a quantum number, and since we study the I z = +1 
sector, isospin mixing can only take place with I > 2 states. In principle there could be mixing 
with rather light / = 2 w + vr states, but as mentioned above this is prevented by neglecting 
disconnected diagrams, which enforce valence flavour structure udss, i.e. 1 = 1 and I z = +1. 

To summarise, for the ao(980) sector I(J P ) = 1(0 + ) mixing due to twisted mass symmetry 
breaking is not expected to cause any problems. This is confirmed by our numerical results (no 
additional unexpected states are observed, the effective mass plateaux quality is good and does 
not seem to be contaminated by mixing; cf. section 

3.1.2 Different twisted mass realizations of the s quark 

In our mixed action setup the s quark can be realized with either a twisted mass term +ifi s 'f5 
or — denoted by s + and s~ , respectively (cf. also section l2TTj) . Consequently, the ss pair 
appearing in our creation operators can be s + s + , s~s~, s + s~ or s~s + . In the continuum limit 
all four choices yield identical results. At finite lattice spacing, however, results are different due 
to discretisation errors. 

For mesons, it is known that using a quark and an antiquark with different twisted mass signs 
significantly reduces discretisation errors [40} |4"T| 142] . With this in mind s + s~ should be the 
optimal choice for the operators ©, Q and (fTUj) . 

It is not clear, whether this mixed realization is optimal also for diquarks (operator (JSJ)). For this 
reason, we also used s + s + (or s~s~ , which yields exactly the same result). Another advantage is 
the possibility to also compute disconnected diagrams (which we plan to do in the near future), 
which is not possible, when using s + s~. 

Performing computations both with s + s + as well as with s + s~ is not only a valuable cross check 
of numerical results, but also provides a first estimate of the magnitude of discretisation errors 
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at our current value of the lattice spacing. In section [J] numerical results are presented and 
discussed in this context. 



3.2 Creation operators, k sector (quantum numbers I(J P ) = l/2(0 + )) 
The expected low-lying spectrum in the n sector (~ 700 MeV) is the following: 

• A two-particle K + ir state. 

- In nature: mass m(K + vr) « 496 MeV + 140 MeV = 636 MeV [I2j . 

— In our lattice setup: masses m{K + tt) ~ m(K) + m(ir) depend on the gauge link 
ensemble and can be read off from Table [2J 

• Possibly a bound k state (might be of quark-antiquark, molecule or diquark-antidiquark 
type), mass m = 682 ± 29 MeV [12]. Such a state has been observed in the lattice study 
reported in [20j, but not in the one in [43]. Note that disconnected contributions are 
neglected in our calculations, like in 



To be able to resolve these low-lying states, we proceed similar as for ao(980) and consider the 
following operators: 

• Molecule type operator: 

qKtt molecule 

= ^2 ( (*( x )75^( x )) (w(x)75«(x)J + (f(x)7 5 d(x)J (d(x)7 5 u(x) 

+ (a(x)75a(x))(a(x)7s«(x))j; (11) 

to be able to check and compare with results of a recent similar lattice tetraquark study 
of k [20], we also consider molecule type operators with 75 replaced by jj and by 77-75. 
Such a structure corresponds to a bound state of a pair of vector mesons (7^) and pair of 
axial vector mesons (7^75), which are significantly heavier than pseudoscalar mesons (75). 
Therefore, we do not expect them to be very helpful to resolve low lying states, which 
is confirmed by our numerical results (cf. section [4731) . They are, however, useful for the 
extraction of excited states. 

• Diquark type operator: 

^diquark = ^ ^ abC S b (^)Cj 5 d C ' T (^ ade d d ' T (^)C^ U e J (12) 

x 

note that 75 diquark flavour combinations [sfZ][uti] and [ss][sit] do not exist, due to the 
Grassmann property of the quark fields, i.e. [uu] = [ss] = 0. Hence, in contrast to the 
molecule operator (llip there is no sum over light quark flavours in the diquark operator 
(fT2|) : as before, we also consider diquark type operators with 75 replaced by 1; since a 
diquark with spin structure 1 is known to be heavier than a diquark with spin structure 
75 [SHI [39], we mainly use it to resolve excited states. 
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Both our numerical results for ao(980) and the above mentioned lattice study of k |20| indicate 
that two-particle scattering states can be resolved with four-quark operators, where all quarks 
are located at the same point (e.g. operators (fTTj) and (fT2"j) ). Therefore, in contrast to our study 
of ao(980) we do not consider operators of two-particle type. This allows to save a significant 
amount of computer time, because two-particle operators require different inversions of the Dirac 
operator (timeslice propagators instead of point propagators; cf. section [3. 3.2[) . 

3.2.1 Mixing due to twisted mass symmetry breaking 

In contrast to the ao(980) sector, mixing introduces additional low- lying states in the k sector, 
which have to be understood and resolved. These additional states have their origin in two- 
particle K + -/r states with / = 3/2 (an 1 = 1/2 kaon and an / = 1 pion can either form an 
1= 1/2 or / = 3/2 state). 

In QCD, where isospin is conserved, these states are degenerate. One can linearly combine 
the degenerate I z € {—1/2, +1/2} kaons and I z = {— 1,0, +1} pions with appropriate Clebsch 
Gordan coefficients to form states with defined isospin 1 = 1/2 and / = 3/2. Note, however, that 
any other linear combination is still an eigenstate of the QCD Hamiltonian, i.e. when discussing 
eigenstates of the Hamiltonian defined isospin is an option, but not a necessity. 

In twisted mass lattice QCD isospin is broken by 0(a), i.e. u and d as well as s + and s~ quarks 
are treated differently. The important consequence in the context of this discussion is that the 
I z = —1/2 kaon s + d is lighter than the I z = +1/2 kaon s + u. Similarly there is a splitting of 
pion masses, where the I z = ±1 pions (du and ud) are lighter than their I z = counterparts 
(uu — dd). While in QCD any linear combination of these kaons and pions is an eigenstate of 
the Hamiltonian, this splitting determines specific linear combinations, which are eigenstates in 
twisted mass lattice QCD: there is a (K+ir) = (s + u+(uu — dd)) state and a (K+ir) = (s + d+du) 
state; the two mesons in the first state are heavier than the two mesons in the second state (cf. 
Table [2]). Note that both combinations have / = 1/2 and / = 3/2 contributions of the same 
order of magnitude, i.e. are not even close to isospin eigenstates. Thus, when determining the 
low lying spectrum, one needs to resolve / = 1/2 as well as / = 3/2 K + ir states. 

To summarise, for the « sector I(J P ) = 1/2(0+) mixing due to twisted mass symmetry breaking 
will double the number of two-particle K + ir states contained in our correlation matrices. This 
theoretical expectation is confirmed by our numerical results (cf. section [5]). 

3.2.2 Different twisted mass realizations of the s quark 

As mentioned in the previous section, we realize the s quark via s + . 

Using s~ would yield on a quantitative level slightly different numerical results. The reason is 
that one would observe a (K + ir) = (s~u + (uu — d)d) state and a (K + ir) = (s~d + du) state, 
i.e. each of the two states contains one "heavy version" of a meson and one "light version" of a 
meson. Of course, in the continuum limit s + and s~ yield identical results. 

The results presented in this paper exclusively correspond to s + . 
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3.3 Computation of correlation matrices 



We compute separate correlation matrices for ao(980) and k, 



c jk (t 2 -h) = ((Oite^cMti)) 



(13) 



where Oj and Ok denote the creation operators ([6]) to (|12p . Technical details of these computa- 
tions are explained in the following. 

3.3.1 Smearing techniques 

To improve the overlap to the low- lying states of interest, we use Gaussian smearing of quark 
fields [444 145j with APE smeared spatial links [16] . Detailed equations can be found e.g. in |47] . 
We use APE smearing parameters «ape — 0.5 and -/Vape = 20. Gaussian smearing is done with 
^Gauss = 0.5 and -/VQauss = 50 for most ensembles. Only for A40.20 we used A^Gauss = 30 instead 
of N G 

auss — 50. For lattice spacing a ~ 0.086 fm these parameters are inside a region, in which 
the overlap between mesonic trial states and the K and D meson is rather large |36j . 

3.3.2 Propagator computation 

For correlation matrix elements (|13p . where both Oj and Ok are molecule and/or diquark op- 
erators, we use point source inversions, i.e. twelve inversions per gauge link configuration and 
quark flavour. This yields point-to-all propagators, which are exact, but which do not allow 
to exploit spatial translational invariance of the correlation matrix elements, to increase their 
statistical precision. In order to reduce correlations, however, we have chosen a random position 
for the source vector for each gauge configuration. 

For correlation matrix elements, where at least one of the operators Oj and Ok is a two-particle 
operator, the situation is different: here the standard one-end trick |48] can be applied twice, 
allowing a stochastic estimation of timeslice-to-all propagators. For each application of the one- 
end-trick we generated an independent stochastic timeslice source with Z 2 x Z 2 noise, where the 
source time slice has been chosen randomly for each gauge configuration. Computing correlation 
matrix elements with stochastic timeslice-to-all propagators is rather efficient, because they 
allow to exploit spatial translational invariance, which in turn reduces gauge noise significantly. 
Moreover, correlations between two two-particle operators require timeslice-to-all propagators, 
which are prohibitively expensive to compute using point source inversions. 

3.4 Analysis of correlation matrices 

To extract energy levels from N x N correlation matrices, we solve the generalised eigenvalue 
problem 



(cf. e.g. [39] and references therein). For a lattice with infinite temporal extension T the eigen- 
values X n (t,to) are proportional to e~ Ent for sufficiently large t, where E n are the energies of 
the N lowest energy eigenstates. 





ra = 0, ...,N- 1 



(14) 
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However, for lattices with periodic finite temporal extensions and sectors, where light two- 
particle states exist (in our case two pseudoscalar mesons; cf. sections [3 . 1 1 and [33]) . the interpre- 
tation of the eigenvalues X n (t,to) is no longer simple. For example a diagonal correlator Cjj(t), 
which is dominated by a two-particle state with energy E n shows the following behaviour for 
< t < T HHEUCEB]: 



ATI 



-E n t 



+ e 



-E n {T-t) 



+ 



e -m 1 t e -m 2 (T-t) +e - 



■m 2 t- mi {T-t) 



(15) 



where m\ and m,2 denote the masses of the corresponding single-particle states and A™ and B™ 
are operator dependent and problem specific constants. The "BJ term" corresponds to the "mi 
particle" traveling forward in time, while the "m2 particle" is traveling backwards in time, and 
vice versa. This term leads to a drastic and characteristic deviation of effective masses from 
their plateaux values at larger temporal separations. For example in Figure [2b this effect is 
clearly visible for t/a> 15. 

In [20] eq. (fT5j) was fitted to the eigenvalues A n (i,io) to extract the energy levels E n . In this 
project, however, several two-particle states with rather different single-particle masses con- 
tribute: in the ao(980) sector K + K and r] s + ir are the relevant states (cf. section [3TT]) ; for k, 
due to twisted mass symmetry breaking, light and heavy kaons and pions need to be considered 
(cf. section I3.2[) . Since a proper treatment of all these two-particle states yields an equation 
significantly more complicated than (115j) with too many parameters to perform stable fits, we 
follow a different route. 



Since A 7 - and B™ are of 



First note that the £?" term in (fT5|) is suppressed by « e - min ( m i> m 2)T_ 
the same order of magnitude, the BJ term becomes irrelevant for sufficiently small t or T — t. 
Hence, we extract the energy levels considering small temporal separations only. We restrict 
all our effective mass analyses to t,T — t <T/4, which seems to be a rather conservative choice. 
Possibly present excited state contributions are taken into account by fitting two exponentials 
to each of the eigenvalues A n (i,io) of interest, fitting range i m i n < t < i max . to, i m i n and t max 
have been chosen such that x 2 /dof<l. Moreover, we varied the values of to, t m in and i max to 
check and confirm the stability of our fitting results. 



4 Numerical results and their interpretation 

4.1 a (980): tetraquark and two-particle operators, ensemble A40.20 

We start by discussing ao(980) (I(J P ) = l/2(0 + )) results obtained using ensemble A40.20 (cf. 
Table [1]). This ensemble with rather small spatial extension (L ~ 1.72 fm) is particularly suited 
to distinguish two-particle states with relative momentum from states with two particles at rest 
and from possibly existing ao(980) tetraquark states (two-particle states with relative momentum 
have a rather large energy because one quantum of momentum p m i n = 2tt/L ~ 720 MeV). 

Figure [2k. shows effective mass plots from a 2 x 2 correlation matrix with a KK molecule 
operator ([6]) and a diquark-antidiquark operator (jSJ), flavour combination s + s~ (cf. section l3.1.2p . 
The corresponding energies extracted from the two plateaus are given in Table H] and they are 
consistent both with the expectation for possibly existing ao(980) tetraquark states and with 
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two-particle K + K and 7] s + tt states, where both particles are at rest (m(K + K) ~ 2m(K) 
1198 MeV; m(rj s + tt) » m(r? s ) + m(vr) rj 1115 MeV; cf. Table©. 



i f 



i % 



(a) 



m(K[k]) + m(K[-k]) 
m(»fc[*])+m(»r[-fc]) 



^ 8 it St ^ 

J } * # f 



10 15 



4.5 
4.0 
3.5 
3.0 
| 2.5 
f 2.0 
1.5 
1.0 
0.5 
0.0 



(b) 



i i 1 

- 1 

; j } 

- • ui| 




m{K[k])+m(K[-k]) 

m(»;,[*])+m(»r[-*]) . 

mass from fit 













(c) 



(d) 









two-particle (A' + A") i — i — i - 




two-particle (r/ s + tt) i — * — i 




molecule (KK) \ — ■ — i 
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'["] 



t[a] 



Figure 2: ao(980) sector, A40.20 ensemble, (a) Effective masses as a function of the temporal 
separation, 2x2 correlation matrix (local operators: KK molecule, diquark-antidiquark, eqs. 
([6]) and (HI)). Horizontal lines indicate the expected two-particle K + K and t] s + tt energy levels, 
(b) 4x4 correlation matrix (local operators: KK molecule, diquark-antidiquark, two-particle 
K + K, two-particle i] s + it, eqs. © to (fTUj) ). (c), (d) Squared eigenvector components of the 
two low-lying states from (b) as a function of the temporal separation. 

Increasing this correlation matrix to 4 x 4 by adding two-particle K + K and r] s + tt operators 
(eqs. ([!]) and ([TO]) ) yields the effective mass results shown in Figure [2b. Two additional states 
are observed with energies given in Tabled! From this 4x4 analysis we conclude the following: 

• We do not observe a third low-lying state around 1000 MeV, even though we provide 
operators, which are of tetraquark type as well as of two-particle type. This suggests that 
the two low- lying states are the expected two-particle K + K and rj s + n states, while no 
additional stable ao(980) tetraquark state is detected. 

• The effective masses of the two low-lying states are of much better quality in Figure [2b 
than in Figure [2b,. We attribute this to the two-particle K + K and rj s + tt operators, 
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which appear to create larger overlap to those states than the tetraquark operators. This 
in turn confirms the interpretation of the two low- lying states as two-particle states. 

• To investigate the overlap in a more quantitative way, we show the squared eigenvector 
components of the two low-lying states in Figure [2)3 and Figure [2U (cf. |37] for a more 
detailed discussion of such eigenvector components). Clearly, the lowest state is of rj s + -it 
type, whereas the second lowest state is of K + K type. On the other hand, the two 
tetraquark operators are essentially irrelevant for resolving those states, i.e. they do not 
seem to contribute any structure, which is not already present in the two-particle operators. 
These eigenvector plots give additional strong support of the above interpretation of the 
two low lying states as two-particle states. 

• The energy of two-particle excitations with one relative quantum of momentum can be 
estimated by 

m(l + 2,p = Pmin ) « yj m(l) 2 + p 2 min + sj m(2) 2 + p 2 min , Pmin = ^. (16) 

Inserting m(K), m(r] s ) and m(ir) from Table [2J yields m(K + K,p = p ra \ n ) ~ 1873 MeV 
and m(r] s + ir,p = p m i n ) ~ 1853 MeV for the A40.20 ensemble. These numbers are con- 
sistent with the effective mass plateaus of the second and third excitation in Figure [2b. 
Consequently, we also interpret them as two-particle states. 

Figure [2^ and Figure [2b also demonstrate an important technical aspect: two-particle states 
can be resolved by tetraquark operators, i.e. two-particle operators are not necessarily needed, 
to extract the full spectrum. Since we are mainly interested in possibly existing states with a 
strong tetraquark component, we restrict the correlation matrices computed for other ensembles 
to four-quark operators (cf. Table [3]). This allows to save a significant amount of computer 
time, because two-particle operators require different inversions of the Dirac operator (timeslice 
propagators instead of point propagators; cf. section [3.3.2p . 

4.2 a (980): tetraquark operators, many ensembles 

We have analysed the six ensembles listed in Tabled] with respect to ao(980) in a similar way 
as explained in section 14. 11 

As already mentioned above the main difference is that this time we exclusively use tetraquark 
operators ([6]) to ([8]), but no two-particle operators ([9]) and (fTOl) . To be able to resolve more 
than two low-lying states, we supplement the molecule operators and the diquark-antidiquark 
operator by versions, where 75 has been replaced by jj and 7^75 (molecule) and by 1 (diquark- 
antidiquark). More detailed information including e.g. smearing parameters, number of gauge 
link configurations, etc. are collected in Table El 

On a qualitative level our findings agree for all ensembles, i.e. are as reported in the previous 
subsection (effective mass plots are collected in Figure [3]): there are always two low- lying states, 
whose masses are consistent with the expected masses of the two-particle K + K and r] s + tt 
states (cf. Figure [Hand Tabled]); higher excitations (the third, forth, etc. extracted state) are 
in all cases significantly heavier and consistent with two-particle excitations with one relative 
quantum of momentum (cf. eq. (|16p ). 
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Figure 3: ao(980) sector, various ensembles, some of them with twisted mass strange quarks 
others with s + s~ . Effective masses as a function of the temporal separation. Horizontal 
lines indicate the expected two-particle K + K and r? s + ir energy levels. 
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ao (980) sector k sector 
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Figure 4: The two lowest energy levels Eq and E\ obtained by our simulations in the ao(980) 
(s + s + only) and in the k sector (cf. also Table UD as a function of m n +. Additionally, the energy 
levels of the non-interacting states are included as solid lines. 





ensemble 


to 


^min 


^max 


E in MeV 


Ei in MeV 


a (980) 


A40.24 


2 


5 


12 


1199(43) 


1343(24) 


(s+s+) 


A30.32 


3 


7 


16 


1078(129) 


1327(66) 




A80.24 


3 


4 


12 


1321(39) 


1408(31) 


a (980) 


A40.20 


1 


5 


12 


1073(48) 


1195(51) 




A40.32 


3 


7 


16 


1098(77) 


1210(40) 




A50.32 


5 


6 


16 


1130(77) 


1236(48) 


K 


A30.32 


3 


7 


16 


888(30) 


1243(72) 




A40.24 


3 


5 


12 


905(47) 


1316(87) 




A80.24 


3 


5 


12 


1060(50) 


1345(150) 



Table 4: The two lowest energy levels Eq and E\ in the ao(980) and in the k sector (cf. also 
Figure H]). 



To summarise, in the lattice setup and ensembles we are studying there is no indication of any 
additional low-lying tetraquark state. 

4.3 k: tetraquark operators, many ensembles 

The analysis for the k sector (I(J P ) = l/2(0 + )) closely parallels the analysis of the ao(980) 
sector presented above. 

We consider correlation matrices containing a Kn molecule operator (jlip and analogue versions 
with 75 replaced by jj and 7^75 as well as an diquark-antidiquark operator ( (I12p and a similar 
operator with 75 replaced by 1). More detailed information including e.g. smearing parameters, 
number of gauge link configurations, etc. are collected in Table El 

As has been explained in section 13.2.11 in twisted mass lattice QCD isospin / is not a quantum 
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number. Therefore, it is not sufficient to only resolve 1=1/2 two-particle K + tt states. One 
has to take into account also mixing with / = 3/2 two-particle K + tt states, i.e. it is necessary 
to resolve these two types of low-lying two-particle states at the same time. 

Effective mass plots for ensembles A30.32, A40.24 and A80.24 (cf. TableU]) are shown in FigureO 
together with the expected energy levels of two-particle K + tt states (obtained via eq. ()16|) and 
the meson masses collected in Table [2]). While effective mass plateaus are consistent with 
these expected two-particle energy levels, there is no indication of any additional low lying 
state, i.e. of a possibly existing bound k state. While this is suggested by experimental data, 
it contradicts the findings of a similar recent lattice study of n |20| . Currently we have no 
explanation for this qualitative discrepancy of two rather similar lattice computations (same 
operators, no disconnected diagrams, similar quark masses). 
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mass from fit 




A40.24 

m(K\k]) + m'(ir[-*]) (light I =1/2,3/2 comb.) 
m{K[k]) +m(7r[-*:]) T (h|eavy I = 1/2,3/2 comb.) 
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m(K[k\) + m(ir[-*]) (heavy / = 1/2,3/2 comb.) 

mass from fit 



5± 
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Figure 5: k sector, various ensembles. Effective masses as a function of the temporal separation. 
Horizontal lines indicate the expected two-particle K + tt energy levels. 



Results for the two lowest energy levels are collected in Figure S] and Table HI 
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5 Conclusions and outlook 



This work represents a first necessary step in the long term project of studying the scalar mesons 
and their properties on the lattice. The main goal of this work was to develop and test those 
techniques that can be effectively exploited for studying the contribution of four-quark operators 
in mesons, especially in the context of the twisted mass formulation of lattice QCD. 

In particular we computed the low- lying spectrum in the ao(980) and k sectors by employing 
trial states designed to have a substantial overlap with both two-particle and possibly existing 
tetraquark states. With our ensembles we did not see additional states beside those that can 
be identified with the expected two-particle spectrum. In fact for all our ensembles we observed 
two low lying states in correspondence with the K + K and rj s + tt thresholds in the ao(980) 
sector and the K + tt (I = 1/2 and I = 3/2) threshold in the k sector. The next states appear 
roughly consistent with excitations of the first quantum of momentum [2-k/L) on top of those 
thresholds. This is somewhat difficult to reconcile with the additional state found by [20] in the 
k, in spite of the rather similar lattice setups. 

We find that the low lying spectrum has essentially exclusively overlap to two-particle trial 
states. This suggests that the states that we see are, indeed, the expected two-particles states 
at the threshold and not tightly bound states either of molecular type or diquark-antidiquark 
type. 

On the basis of this, we can conclude that either our choice of operators has negligible overlap 
with the wave function of the resonances oo(980) and k, or that our volumes are not large enough 
to identify those states. 

These conclusions can be strengthened by studying more volumes, by introducing twisted bound- 
ary conditions [51] and by studying further trial states of different type. As for the latter, it 
will be crucial to combine four quarks with traditional quark-antiquark operators, but discon- 
nected diagrams will be necessary for that. As for the volume dependence, we plan to use the 
finite volume formulae of Liischer \52\ \T3[ [53"1 [T^l [T5] and their extensions to multiple channels 
developed in [5U E5J ESI E7] . At present, our limited number of volumes is insufficient for such 
an analysis. Corresponding computations are in progress. 

Another possible development consists in studying four-quark states that include the charm 
quark. This is a natural extension thanks also to the presence of a dynamical charm quark 
in the ETMC gauge configurations. This direction is also being explored in particular in the 
context of the tetraquark candidates D* and D*-^ . 
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